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Abstract:  In  this  investigation,  a  two-dimensional  time-domain  closed-form  mathematical  model  for  the  hydrodynamic  forces  on  the 
upstream  vertical  face  of  a  given  rigid  dam  subjected  to  a  specified  horizontal  ground  motion  accelerogram  was  developed.  The  model 
includes  the  absorption  of  energy  at  the  elastic  reservoir  bottom,  characterized  by  the  impedance  ratio  of  the  sub-bottom  materials  with 
respect  to  water  (a).  The  formulated  boundary-value  problem  is  solved  in  Laplace’s  domain  and  subsequently  transformed  back  to  the 
time  domain.  Response  spectra  for  the  hydrodynamic  base  shear  force  and  overturning  moment  are  constructed  for  extreme  values  of  the 
parameter  a.  It  is  found  that,  frequently,  including  the  solid-foundation  elasticity  in  the  reservoir  model  attenuates  the  resultant  hydro- 
dynamic  forces  on  a  rigid  barrier,  as  compared  to  the  results  for  the  case  of  a  rigid  reservoir  foundation.  In  this  case,  the  elasticity  of  the 
sub-bottom  materials  constitutes  an  effective  energy  dissipating  mechanism  (radiation  damping).  By  contrast,  for  sub-bottom  materials 
with  less-than-water  impedance,  amplification  of  the  effective  earthquake  forces  is  obtained,  as  compared  to  the  results  for  the  case  of  a 
rigid  reservoir  foundation. 
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Introduction 

Realistic  assessment  of  the  response  of  dam-reservoir-foundation 
systems  to  earthquake  excitation  requires  appropriate  modeling  of 
the  sub-bottom  absorption  of  energy.  The  acoustic  impedance  (or, 
equivalently,  the  coefficient  of  reflection)  of  the  sub-bottom  ma¬ 
terials  in  relation  to  water  has  been  found  to  adequately  represent 
the  corresponding  equivalent  damping  in  the  boundary  conditions 
of  the  governing  equations  of  motion  for  the  water  medium  in  the 
frequency  domain  (Hall  and  Chopra  1982;  Fenves  and  Chopra 
1985b).  However,  the  exact  interaction  of  the  reservoir  with  its 
foundation  remains  difficult  to  characterize. 

When  a  water  pressure  wave  traveling  vertically  downward 
strikes  the  bottom  boundary  surface,  part  of  the  propagating  en¬ 
ergy  is  reflected  back  into  the  reservoir  and  another  part  is  trans¬ 
mitted  to  the  foundation  medium  on  account  of  its  finite  stiffness. 
Frequently,  this  refracted  energy  may  be  assumed  to  continue 
propagating  vertically  downward  toward  infinity,  actually  consti¬ 
tuting  an  equivalent  energy  dissipating  mechanism  introducing 
radiation  damping  that  effectively  reduces  the  response  of  the 
system.  On  the  other  hand,  if  the  sub-bottom  materials  contain 
large  amounts  of  constituents  with  lower-than-water  impedance,  a 
substantial  part  of  the  incident  pressure-wave  energy  may  be  re¬ 
flected  back  into  the  reservoir  as  rarefaction  waves,  which  in  turn 
become  reflected  compression  waves  at  the  upper  boundary  sur- 
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face  (on  air),  which  may  reinforce  the  continuous  input  of  pres¬ 
sure  from  the  exciting  dam  motion.  In  other  words,  it  is  possible 
that  the  interaction  of  the  reservoir  with  the  sub-bottom  materials 
effectively  increases  the  magnitude  of  the  system  response 
(Cheng  1986;  Bougacha  and  Tassoulas  1991).  The  resulting  ef¬ 
fect,  as  compared  to  the  response  for  a  rigid  foundation,  depends 
on  both  the  actual  impedance  ratio  of  the  foundation  with  respect 
to  water  and  on  the  specific  characteristics  of  the  input  ground 
motion  acceleration,  particularly  on  the  content  of  large  pulses. 

This  investigation  presents  a  mathematical  model  to  evaluate 
the  force  effect  by  the  impounded  water  on  retaining  concrete 
gravity  dams  upon  earthquake  excitation,  in  the  time  domain, 
including  the  interaction  of  the  reservoir  with  the  sub-bottom  ma¬ 
terials.  The  time  history  of  the  hydrodynamic  pressure  field  on  the 
upstream  face  of  the  dam  and  the  associated  resultant  lateral  force 
and  overturning  moment  are  evaluated,  allowing  the  construction 
of  the  corresponding  response  spectra  for  dam-reservoir-sub- 
bottom  systems. 

The  analytical  model  is  two  dimensional  (which  is  very  appro¬ 
priate  for  straight  gravity  dams)  and  the  dam  is  considered  as  a 
rigid  barrier.  The  reservoir  is  assumed  infinitely  long  in  the  up¬ 
stream  direction.  Water  is  assumed  a  homogeneous,  inviscid,  and 
linearly  compressible  fluid,  and  its  flow  is  taken  as  irrotational. 
The  boundary  condition  at  the  horizontal  bottom  of  the  reservoir 
is  formulated  accounting  for  the  elasticity  of  the  underlying  me¬ 
dium.  The  sub-bottom  is  represented  by  means  of  a  characteristic 
value  of  its  impedance  ratio  with  respect  to  water.  The  exciting 
input  ground  motion  acceleration  is  assumed  to  be  acting  in  the 
horizontal  direction  (parallel  to  the  reservoir  bottom)  only  and 
the  initial  conditions  of  the  system  are  assumed  to  be  at  rest. 
Hydrostatic  forces  and  the  secondary  effect  of  gravity  waves  are 
ignored. 

The  direct  analysis  in  the  time  domain  offers  some  definite 
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Fig.  1.  Two-dimensional  analytical  model  for  a  rigid  dam,  the  impounded  water,  and  the  reservoir  bottom 


advantages  over  the  alternative  conventional  techniques  in  the 
frequency  domain.  First,  the  resulting  models  can  easily  be  ex¬ 
tended  to  represent  the  behavior  of  flexible  dams  exhibiting  ma¬ 
terial  nonlinearities.  Nonlinear  frequency-domain  analysis,  even 
when  possible  at  all,  is  extremely  cumbersome.  Second,  analysis 
in  the  time  domain  is  agreeable  to  the  intuition  of  engineers.  In 
fact,  direct  time-domain  analysis  using  a  closed-form  solution  al¬ 
lows  practical  and  efficient  assessment  of  the  hydrodynamic 
forces  on  barriers  in  reservoir-dam  systems  subjected  to  general 
earthquake  excitations  (much  of  the  seismic  frequency-domain 
researches  in  dam  engineering  to  date  have  been  limited  to  har¬ 
monic  earthquakes  (Bouaanani  et  al.  2003;  Fenves  and  Chopra 
1985a;  Gogoi  and  Maity  2006).  The  corresponding  numerical 
analyses  in  the  frequency  domain  are  more  intricate,  involving  the 
use  of  fast-Fourier-transform  algorithms  (Fenves  and  Chopra 
1983,  1984a, b).  Finally,  but  most  important,  the  simple  closed- 
form  solution  in  the  time  domain  allows  expedient  and  realistic 
assessment  of  the  frequent  economical  advantages  of  including 
bottom  absorption  in  the  engineering  estimation  of  earthquake 
hydrodynamic  forces  on  dams,  including  the  efficient  generation 
of  response  spectra,  for  a  given  site  and  design  ground  motion. 
Previous  investigations  in  the  time  domain  invariably  implement 
numerical  methods  of  solution  (Weber  1994). 


Mathematical  Models 
Problem  Formulation 

The  physical  behavior  of  the  dam-reservoir-sub-bottom  system  is 
described  mathematically  by  the  two-dimensional  wave  equation 
in  terms  of  the  velocity  potential  for  the  water  particles  in  the 
reservoir,  c|>  (Kotsubo  1959) 

=  °  (1) 

where  V2=Laplacian  operator  d1  /  dx2 + d1  /  dz2  and  c=velocity  of 
propagation  of  the  pressure  waves  in  the  water. 

Fig.  1  shows  schematically  the  boundary  conditions  at  the  four 
edges  of  the  water  domain,  identified  as  (i)-(iv).  Boundary  con¬ 
dition  (i)  refers  to  the  bottom  of  the  reservoir  and  represents  the 
absence  of  the  vertical  component  of  acceleration  for  the  water 


particles  in  contact  with  the  solid  bottom  material  since  the  input 
ground  motion  is  assumed  to  act  exclusively  in  the  horizontal 
direction. 

In  the  frequency  domain,  the  gradient  of  the  Fourier  transform 
of  the  pressure  field  in  the  direction  of  the  local  normal  to  the 
contact  surface  at  the  bottom  consists  of  two  distinct  terms.  One 
term  due  to  the  free-field  ground  acceleration  at  the  site,  with 
local  normal  component  an,  and  another  term  due  to  the  static 
interaction  between  the  solid  and  the  fluid  media  upon  reflection 
and  refraction  of  the  pressure  waves  at  the  interface  between 
these  two  phases  owing  to  their  different  elastic  properties 
(Rosenblueth  1968;  Hall  and  Chopra  1982).  The  corresponding 
relationship  may  be  expressed  mathematically  as 

<?P  (w) 

— - =  -  p  •  a„  +  ( iw )  ■  q  ■  P(w)  (2a) 

dn 

where  P(w)= Fourier  transform  of  the  pressure  field  in  the  water; 
p=density  of  water;  w= frequency-domain  variable;  i=  imaginary 
unit;  and  q  is  given  by 


q= -  (2b) 

Pr  '  Cr 

in  which  p,.= density  of  the  sub-bottom  material  and  cr= velocity 
of  propagation  of  the  pressure  waves  through  the  same  medium, 
i.e.,  the  product  (q-c)  is  the  impedance  ratio  of  water  with  respect 
to  the  foundation  material.  Alternately,  the  product  (q  ■  c)  may  also 
be  expressed  in  terms  of  the  coefficient  of  reflection  of  the  sub¬ 
bottom  material  with  respect  to  water  ((3,.),  defined  for  the  pure- 
interaction  problem  (i.e.,  when  there  is  no  free-field  ground 
motion  acceleration),  as 


I  -  Pr  1 

q  ■  C  = - =  — 

1  +  Pr  a 


(2c) 


where  a = impedance  ratio  of  the  sub-bottom  material  with  respect 
to  water,  i.e. 


P  r'Cr 
a  = - 

P  '  c 


(2d) 


Since  the  vertical  component  of  ground  motion  acceleration  (an) 
is  zero  in  the  problem  being  formulated,  boundary  condition  (i), 
expressed  by  Eq.  (2a)  in  the  frequency  domain,  reduces  to 
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(9P(w) 

dn 


=  ( iw )  ■  q  ■  P(vv) 


which  may  be  written  in  the  time  domain  as 


dp  dp  dp  dp 

—  -  q  ■  —  =  —  -  q  ■  —  =  0,  at  z  =  0 
dn  dt  dz  dt 


(3a) 


(3b) 


where  p(x,z,t)  = hydrodynamic  pressure  field,  given  by  (New- 
mark  and  Rosenblueth  1971) 


p(x,z,t)  =  p 


(?4> 

dt 


(4) 


Therefore,  boundary  condition  (i)  in  the  time  domain  becomes 


dz 


=  q- 


(9<f> 

dt 


(5) 


Eq.  (5)  is  the  distinctive  boundary  condition  for  the  time  domain 
analysis  of  a  dam-reservoir-flexible-subbottom  system  subjected 
to  earthquake  excitation.  The  other  three  boundary  conditions  re¬ 
main  the  same  as  for  the  associated  rigid  foundation  problem 
(Kotsubo  1959),  i.e., 


(ii) 


The  horizontal  velocity  of  the  water  particles  in  direct  con¬ 
tact  with  the  upstream  face  of  the  dam  is  given  by  the 
function  for  the  ground  motion  velocity.  This  is  expressed 
mathematically  as 


(9<j> 

dx 


=  ve(t) 


(6) 


(iii) 


The  hydrodynamic  pressure  at  the  upper  edge  of  the  water 
domain  is  zero  (ignoring  the  effect  of  gravity  waves).  Or, 
expressed  mathematically,  in  terms  of  the  water  velocity 
potential,  one  has 


d§ 

dt 


=  0 


(V) 


(iv) 


Very  far  away  from  the  direct  disturbance  of  the  water 
medium,  all  hydrodynamic  responses  of  the  water  should 
tend  asymptotically  to  zero  (radiation  condition).  This  is 
expressed  mathematically  as 


lim  cf>  =  0 


(8) 


To  complete  the  formulation  of  the  problem  in  the  time  domain, 
in  addition  to  these  boundary  conditions,  zero  initial  conditions 
are  assumed  for  the  system,  i.e.,  in  terms  of  the  water  velocity 
potential,  one  has 


4>(0)  =  const  (=0) 


and 


di |> 
dt 


(0)  =  0 


(9a) 


(9b) 


Laplace  Transformation  of  the  Problem 

Taking  the  Laplace  transformation  of  the  wave  Eq.  (1)  with  re¬ 
spect  to  time,  the  governing  differential  equation  becomes  the 
following  Helmholtz  equation 


V2<J>  -  —  •  $  =  0 


(10) 


where  <l>=T(tp) = Laplace  transformation  of  the  water  velocity  po¬ 
tential,  and  s  is  the  basic  variable  in  Laplace’s  domain.  Boundary 
conditions  (i)-(iv)  may  also  be  correspondingly  transformed  to 
give 


dz 


=  q  ■  s  ■  (4>)|.=0 


z=0 


d<S> 

dx 


=m«] 


and 


=  o 


lim  <J>  =  0 

X — tco 


(11a) 

(11  b) 

(11c) 

(ii  d) 


Problem  Transformation  within  Laplace’s  Domain 

Looking  for  homogeneous  boundary  conditions  with  respect  to 
the  independent  variable  z,  let  the  transform  of  the  unknown  ve¬ 
locity  potential  be  defined  in  terms  of  a  new  field  in  Laplace’s 
domain,  'I7,  according  to 


<j)  =  e -q  s-th-z) .  ^ 


(12) 


When  Eq.  (12)  is  substituted  in  Eq.  (10),  the  new  governing  dif¬ 
ferential  equation  is  obtained  as 


V2T7  +  2-q-s- - +  (q  ■  s)2  ■  I  1  - 

dz 


1 


(q-c)2 


■¥  =  0  (13) 


and  the  boundary  conditions  expressed  by  Eqs.  (11a)  and  (11c) 
become 


dV 

dz 


=  0 


z=0 


and 


respectively. 


*U=o 


(14a) 


(14  b) 


Problem  Solution  in  Laplace’s  Domain 

This  boundary-value  problem  is  conveniently  solved  for  the  field 
1I7  in  Laplace’s  s-domain  by  separation  of  variables.  By  assuming 
1P=V(jc,5) -Z(z,s),  Eq.  (13)  may  be  written  as  (Wylie  1975) 


X"  s2\  Z"  Z'  ,  ,, 

tr?M7+2'‘,"'7+‘' ,s  1=0 


(15) 


where  Z'  =  dZI dz,  Z" =d1ZI dz2,  and  X'^cfiXI dx2.  One  way  to  sat¬ 
isfy  Eq.  (15)  is  by  having  the  first  parenthesis  in  this  equation 
equal  to  a  positive  constant  X2,  while  the  second  parenthesis  in 
the  equation  equal  to  the  negative  of  the  same  constant,  i.e. 


—  -  —  =  X2 
V  c2 


(16a) 


and 
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eta 


-f\m  =  Ki-h^O--m+l)--  (m  53  0)  (21) 

independently  of  s.  It  will  be  subsequently  shown  that  the  contri¬ 
butions  to  the  solution  from  very  large  values  of  s  are  negligible. 
Therefore,  the  eigenfunctions  generated  by  Eq.  (17)  may  be  writ¬ 
ten  as 

Zm(z,s)  =  Am(s)  ■  e~qsz  ■  |cos(X,„  ■  z)  +  •  sin(Xm  ■  z)j 

(22) 

where  Xm  is  given  by  Eq.  (21). 

Now,  Eq.  (16a)  renders  the  following  ordinary  differential 
equation: 


Fig.  2.  Graphical  representation  of  the  numerical  solution  for  the 
eigenproblem  formulated  in  Eq.  (20) 


—  +  2  ■  q  ■  s  ■  —  +  q2  ■  s2  =  -\2  (1 6b) 

The  general  solution  of  Eq.  (16h)  is 

Z(z,  s)  =  e~q  s  z  ■  [A(s)  •  cos(X  •  z)  +  B(s)  ■  sin(X  ■  z)]  (17) 

Imposing  boundary  condition  (14a)  into  Eq.  (17)  leads  to  the 
relation 

-  q  ■  s  ■  A(s)  +  X  ■  B(s)  =  0  (18a) 

while  imposing  boundary  condition  (14 b)  into  Eq.  (17)  leads  to 
the  relation 


A( s)  •  cos(X  ■  h)  +  B(s)  ■  sin(X  ■  h)  =  0  ( 1 8h) 

Putting  Eqs.  (18a)  and  (18h)  into  a  single  matrix  equation  gives 
the  following  eigenproblem: 

'  X  1  Uwl  |0] 

cos (X  ■  h)  sin(X  •  h)  J  [s(s)J  [oj 

A  necessary  condition  for  this  algebraic  system  of  homogeneous 
equations  to  produce  a  nontrivial  solution  is  that  the  determinant 
of  the  coefficient  matrix  be  zero.  Such  condition  leads  to  the 
following  characteristic  equation  (which  in  this  case  is  transcen¬ 
dental): 


tan  t)  —  — - ■  t)  (20) 

q  ■  s  ■  h 

where  "q  =  X  -  /7 .  The  solution  of  Eq.  (20)  yields  an  infinite  number 
of  eigenvalues,  as  indicated  schematically  in  Fig.  2.  The  solid 
curves  in  the  figure  represent  the  left  side  of  Eq.  (20),  whereas  the 
dotted  straight  line  represents  the  right  side  of  the  same  equation. 
The  abscissas  of  the  points  of  intersection  of  the  straight  line  with 
the  set  of  curves  (highlighted  in  the  figure)  are  the  solutions  to 
Eq.  (20),  which  must  be  determined  numerically.  In  principle, 
these  eigenvalues  can  be  determined  numerically  with  any  desired 
level  of  accuracy.  However,  it  is  clear  from  the  figure  that,  with 
the  exception  of  the  case  in  which  5  is  “very  large,”  after  a  few 
initial  eigenvalues,  the  mth  root  of  Eq.  (20)  is  given  approxi¬ 
mately  by 


X"  - 


■X=0 


(23) 


with  general  solution 


Xjx,s)  =  Cjs)  ■  +  Dm(s )  •  eW (24) 


However,  the  second  term  in  Eq.  (24)  diverges  for  very  large 
values  of  the  abscissa  x.  Consequently,  on  the  basis  of  the  radia¬ 
tion  boundary  condition  expressed  by  Eq.  (lit/),  we  conclude  that 

Dm(s)= 0. 

Combining  Eqs.  (22)  and  (24)  and  letting  Cm(s)  be  absorbed 
by  Am(s),  one  obtains  the  following  eigenfunction  associated  with 
the  eigenvalue  Xnl : 


' Vm(x,z,s)=Ajs )  ■  e  qs'z  ■  cos(X,„  ■  z )  +  ^  ■  sin(\,„  ■  z) 

\  X,„ 


.  g-VVKP'/cA-*  (25) 

which  is  immediately  transformed  using  the  definition  given  by 
Eq.  (12)  as 


®m(x,z,s)  =  e-q  '  (h-z>  ■  Vm(x,z,s)  (26) 

Finally,  the  general  solution  field  in  Laplace’s  domain  is  the  linear 
combination  of  the  infinite  number  of  eigenfunctions  expressed 
by  Eq.  (26)  to  span  the  entire  space  of  possible  functions 
<f>(x,z,s),  i.e. 


<P(x,z,s)  =  e  q  s  h  ■  2  Ajs)  •  f  cos(\„,  •  z)  +  ^  ■  sin(X„,  ■  z) 

■  g-^n+^lc1)^  (2  7) 

To  proceed  any  further,  one  needs  to  use  the  remaining  boundary 
condition  (11  b)  to  evaluate  the  coefficients  Am(s)  for  a  specific 
input  ground  velocity. 

Unit-Impulse  Response  Function 

The  unit-impulse  response  function  represents  the  hydrodynamic- 
pressure  response  in  the  water  when  the  input  ground  velocity  is 
the  Heaviside  unit-step  function  at  the  initial  time,  or,  equiva¬ 
lently,  when  the  input  ground  acceleration  is  a  Dirac  delta  func¬ 
tion  of  unit  intensity  applied  at  the  origin  of  time.  Imposing 
boundary  condition  (llh)  on  Eq.  (27),  one  obtains 
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dx 


=  -e-^h-JJAm(s)-\K+-2 


q-s 


1 


cos(\m  •  z)  +  —  ■  sin(X„,  •  z)  =  -  -  (28) 


where  the  last  term  in  Eq.  (28)  is  the  Laplace  transform  of  the 
negative  of  the  Heaviside  unit-step  function  at  the  initial  time. 

Any  coefficient  A,-(s)  is  obtained  by  multiplying  each  side  of 
Eq.  (28)  by  cos (X,-z),  integrating  with  respect  to  the  variable  z 
within  the  domain  (-h,+h)  and  considering  the  known  values  of 
the  following  integrals: 


f 

J  -I 


cos (\m  •  z)  •  cos (X,  •  z)dz  =  h  • 


f 


sin(Xm  ■  z)  •  cos(X,-  ■  z)dz  =  0 


(29a) 


(29  b) 


and 


c-rl 


e~(x!c )■  V(X.m  •  c)2+s2- 
_  V(X„,  ■  c)2  +  s2 


=  c  ■  i0[X,„  ■  \j(c  ■  t)2-x2l  c  •  t  >  x 

(33  a) 


The  second  inverse  Laplace  transform  in  Eq.  (32)  may  be  ob¬ 
tained  using  Eq.  (33a)  and  the  fact  that 


dy(t) 

dt 


=  s  ■  7t«p(0]  -  ‘P(0+)  =  *  •  7t<p(f)]  (33 b) 


and  therefore 


rl{s  ■  7t<p(01}  = 


d<p(t) 

dt 


(33  c) 


where  tp(/)  is  an  arbitrary  “well-behaved”  function  of  the  time  t. 

The  complete  Laplace  inversion  of  Eq.  (32)  gives  the  unit- 
impulse  response  function  as  (Bowman  1958) 


f 

J  -i 


cos  (X,-  •  z)dz  =  ■ 


Ah  ■  sin^(2  ■  i  +  1)  ■  — 
(2  •  i  +  1)  ■  TT 


Ah{-  I)1' 

(2  ■  i  +  1)  •  TT 

(29c) 


where  8m,=Kronecker’s  delta.  After  some  algebraic  manipula¬ 
tions,  the  result  is 


Ms)  = 


1  2  sin(X/  •  h > 

'^7 


X,  ■  h  ■  s 


(30) 


Inserting  Eq.  (30),  evaluated  for  i=m ,  into  Eq.  (27),  the  Laplace 
transform  of  the  velocity  potential  is  obtained  as 


F(x,z,t)  =  Apc^j 


(-  D" 


_0  (2  ■  m  +  1)  ■  tt 


jcos(X„,  ■  z)  •  JQ[km  ■  V(c  ■  tY-x2] 

c  ■  t  sin(X„,  •  z)  r  n — ^ 1 

—  ■  ...  ,2  ,  •  ^i[xm  •  V(c  ■  ty  -x2] 

a  V(c  •  t)2-x~ 


c  ■  t>  X 


(34) 


which  simplifies  considerably  when  evaluated  on  the  upstream 
face  of  the  dam  (jc=0)  into 


<t>(x,z,s)  =  2 


1 


2  sin(Xm  ■  h) 


m=  0  '  h  '  S 


X,v  +  ^ 


■  |^cos(X„,  •  z)  +  ■  sin(Xm  •  z)  j  ■  e  '1k+(s~Ic~>* 

(31) 

Notice  that  this  function  goes  to  zero  rapidly  as  s  grows,  confirm¬ 
ing  the  earlier  assertion  that  the  region  of  very  large  values  of  5 
contributes  negligibly  to  the  response.  Using  Eq.  (31),  the  unit- 
impulse  response  function  F(x,z,t)  is  obtained  as 

F(x,z,t )  =  p  ■  7 ~*(j  ■  3>) 

2  sin(Xm  •  h ) 


=  p-2 

m= 0 


X,„  •  h 


!j. 


cos (Xm-z)  ■  r' 


q  ■  sin(xm  ■  z)  j 


■  .  f J-VXm+(s2lc 2)x 


yjX^+is^/c^-x 


K  +  c2 


\ 


Xf,7  +  2 


/ 


(32) 


where  7^*(  )=the  inverse  Laplace  transform  operator. 

The  first  inverse  Laplace  transform  in  Eq.  (32)  is  given  by 
(Zwillinger  1996) 


,  ,  4pCx,  (-1)" 


Tt  «= o  (2  • m  +  !) 


I  c  •  t 
cos(-riffl  ■  0  ■  J0I  T)m  ■  — 


1 


-  -  •  sin(-ri„1  ■  0  ■  7d  T)m  ■  —  f  •  H(t) 


c  ■  t 


(35a) 


where  H(t)= Heaviside  unit-step  function  at  the  origin  of  time; 
i= 0,1,  =Bessel  function  of  the  first  kind  and  ;th  order;  and  £ 
=z//i= normalized  ordinate  of  the  point  under  consideration  along 
the  height  of  the  dam  face.  Function  F°(£ ,  t )  reduces  to  the  solu¬ 
tion  in  the  time  domain  previously  reported  by  Kotsubo  (1959) 
for  a  rigid  foundation  (a— >°°),  i.e. 


,  4pCx,  (-1)" 

nu)=—  2 


V  m= 0  (2  ■'«+!) 


c  •  t 


COs(tlm  ■  0  ■  y0  Tlm  ■  —  ,  f>0 


(35  b) 


Solution  for  Input  Ground  Motion  Accelerogram 


The  seismic  hydrodynamic  pressure  field  on  the  upstream  face  of 
the  dam  when  the  system  is  subjected  to  ground  motion  accelera¬ 
tion  is  readily  computed  from  the  convolution  of  the  unit-impulse 
response  function  with  the  time  history  of  the  ground  motion 
acceleration,  a?(f)  (Clough  and  Penzien  1993) 
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U)=  f 

J  0 


f’ 

Jo 


P°(&)=  I  ag(t  -  T)  ■  F°(i„Odr  =  y  ■  g  ■  |  T(t-j)  ■  F°(£„j)dj 

(36) 

where  p°(£,f)=hydrodynamic  pressure  on  the  upstream  face  of 
the  dam  as  a  function  of  time;  T =nondimensional  input  accelera¬ 
tion  function;  and  y  ■  g  =  peak  ground  motion  acceleration,  as  a 
fraction  (y)  of  gravity  (g). 

Inserting  Eq.  (35a)  into  Eq.  (36),  one  obtains  the 
hydrodynamic-pressure  response  as 


P°(.i,t)  =  4 


■  c  ■  y 


TT 


■2 


(-  D" 


=0  (2  ■  m  +  1) 


cos(ti,„  •  0  ■  A0m(t )  -  -  •  sin(Ti,„  •  Q  ■  Alm(t) 
a 


where  w0=unit  weight  of  water  and 


Aim{t)  =  J  r(t  -  t)  •  ‘  =  0,1 


(37  a) 


(37  b) 


Alternatively,  Eq.  (37a)  may  be  normalized,  for  a  nondimensional 
solution  for  the  hydrodynamic  pressure,  as 


T  °(i,t)  = 


=  -•2 


(-  ir 


Wo-c-y  -it  m=0  (2-  m  +  1) 


cos(nqm  •  0  ■  A0m(t )  -  -  •  sin(in,„  •  0  ■ Alm(t ) 
a 


(38) 


Resultant  Hydrodynamic  Force 


The  resultant  earthquake  hydrodynamic  force  on  the  upstream 
face  of  the  dam  P°(t)  is  obtained  by  integrating  the  pressure  held 
along  the  height  of  the  reservoir,  i.e. 


P°(t)  =  h 


■  [  P°(C,t)dt 

Jo 


Inserting  Eq.  (37)  into  Eq.  (39)  and  recognizing  that 


f 

J  o 


cos('T)m  •  Qfite  = 


(-  ir 


/  1\ 

\  +  2/ 

•  77 

(39) 


(40a) 


and 


_Jr(f)__  32  v 
fin  1  2  7  f  ■'  • 


,=o  (2  •  m  +  1) 


a 


(41*) 


where  f=c/(4-h)  is  the  fundamental  frequency  of  natural  vibra¬ 
tions  of  the  reservoir  (Newmark  and  Rosenblueth  1971). 

Resultant  Hydrodynamic  Overturning  Moment 

The  resultant  earthquake  hydrodynamic  overturning  moment  on 
the  base  of  the  dam  M°(t)  is  obtained  by  integrating  along  the 
height  of  the  reservoir  the  forces  exerted  by  the  pressure  held 
factored  by  their  corresponding  lever  arms  with  respect  to  the 
base,  i.e. 


M°(t)  =  h2  ■  £-p°(l,t)dt 


Inserting  Eq.  (37)  into  Eq.  (42)  and  recognizing  that 


■f 

Jo 


(42) 


jv. 

j  0 


■  cos(ti,„  ■  0 dq  =  — - — - 

/  o  (2  ■  m  +  1 )  ■  tt 


(-i  r- 


(2  ■  m  +  1)  •  ir/ 

(43a) 


and 


r  e  si 

Jo 


si n(Vm-C)dC  =  (~  !)" 


(2  •  m  +  1)  ■  it 


(43*) 


the  resultant  hydrodynamic  overturning  moment  on  the  dam  base 
is  obtained  as 


. .  wn  ■  c  ■  y  •  h  ^ 
M°(t)  =  16  •  — - y - •  2 


(-  l)m 

"  m=0  (2  ’  m  +  l)3 

-  A0m(t)  •  (1  -  (-  1)'"  ■  O  -  -  ■  (-  1  r-Alm(t) 

OL 

(44  a) 

or,  normalizing,  for  a  nondimensional  solution,  Eq.  (44a)  may  be 
written  as 


M°(t)  64  v 

>  3  3  ■  y  ■  f ' 

Wn  ■  n  TT 


1 


=o(2- 


l)3 


^0 m(t)  '  [(“  1)HI+I  +  Tim]  -  ■■AlmW 
a 


(44  b) 


1 


[  sin(T)m  ■  Qdi  «  — 

J0  Pn 


the  resultant  hydrodynamic  force  is  obtained  as 

oo 

w„  -  c  ■  y  ■  h  ^  1 


P°(t)  =  8  •  ■ 


-■2 


i=o  (2  ■  m  +  1  )2 


1 


A0m(t)  ~  0m  ■  Alm(t) 
a 


(40*) 


(41a) 


or,  normalizing,  for  a  nondimensional  solution,  Eq.  (41a)  may  be 
written  as 


Case  Study 

A  simple  case  study  provides  insight  into  the  influence  of  the 
foundation  impedance  ratio  with  respect  to  water  (a)  on  the  re¬ 
sponse  of  rigid  dams  to  horizontal  ground  motion  acceleration. 
During  this  study,  three  extreme  values  of  a  (0.1,  1.0,  le+6)  were 
considered  to  emphasize  the  effect  of  including  the  elasticity  of 
the  sub-bottom  materials  in  the  model.  The  hrst  value  (a=0.1) 
represents  a  “highly  gaseous”  sub-bottom.  These  are  hypothetical 
materials  with  a  large  volumetric  content  of  undissolved  gases  or 
any  other  components  that  constitute  a  phase  with  substantial 
less-than- water  impedance.  The  second  value  (a=l)  represents  a 
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Fig.  3.  Unit-impulse  response  function  for  the  reservoir  hydrody¬ 
namic  pressure  on  the  vertical  face  of  a  rigid  dam  facing  a  headwater 
height  of  h  =  65  m,  at  a  nondimensional  height  of  £  =  0.40,  for  several 
values  of  the  bottom  impedance  ratio  a 


type  of  reservoir  bottom  that  is  just  an  extension  of  the  water 
medium  down  toward  infinity.  Liquefied  reservoir  bottoms  with 
solid  particles  in  suspension  would  have  a  values  in  the  upper 
vicinity  of  a=  1.  The  last  value  (a  =  le+6)  represents  a  very  rigid 
sub-bottom,  as  in  models  that  do  not  include  interaction  with  the 
sub-bottom  materials. 

Fig.  3  shows  the  normalized  unit-impulse  response  function 
for  the  reservoir  hydrodynamic  pressure  on  the  vertical  face  of  a 
rigid  dam  facing  a  headwater  height  of  h  =  65  m,  at  a  nondimen¬ 
sional  height  of  £=0.40,  for  the  three  extreme  values  of  a.  The 
function  F°(£,f)/(p-c)  is  given  by  Eq.  (35 a).  Notice  the  effect  of 
including  the  parameter  a  in  the  formulation.  At  this  scale,  the 
plot  does  not  show  much  difference  between  the  functions  for  a 
rigid  and  for  elastic  solid  sub-bottoms.  But  the  function  is  dra¬ 
matically  magnified  for  highly  gaseous  sub-bottom  materials. 
This  extreme  contrast  is  indicative  of  hydrodynamic-force  ampli¬ 
fications  when  a  is  smaller  than  1  (more  and  more  so  as  a  de¬ 
creases  toward  its  lower  bound,  zero),  particularly,  if  the  input 
ground  motion  accelerogram  contains  large  pulses. 

Fig.  4  shows  the  input  accelerogram  selected  in  the  study.  This 
is  a  short-duration  (slightly  less  than  2  s)  segment  of  the  recorded 
and  base-corrected  N21E  component  of  the  Taft  earthquake  (July 
1952).  The  peak  ground  acceleration  is  only  0. 1 6g,  but  the  record 
contains  a  large  pulse  of  input  energy. 

Fig.  5(a)  shows  the  normalized  base  shear-force  response  of  a 
rigid  dam  facing  a  headwater  height  of  h  =  65  m,  when  subjected 
to  the  selected  input  earthquake,  for  the  three  extreme  values  of 
a.  The  function  110  is  given  by  Eq.  (41  b).  As  expected,  the  hy¬ 
drodynamic  force  response  is  dramatically  magnified  for  the 
highly  gaseous  sub-bottom  materials. 

To  appreciate  the  effect  of  including  the  a  parameter  for  solid- 
elastic  sub-bottoms  in  the  reservoir  model,  the  functions  for 
a  =  le  +  6  (rigid  sub-bottom  materials)  and  a=l  (water  extension) 
are  plotted  again  (at  the  appropriate  scale)  in  Fig.  5(b).  Notice 
that,  generally,  there  is  an  attenuation  of  force  response  when  the 
elasticity  of  the  reservoir  foundation  is  included  in  the  model.  The 
parameter  a  clearly  provides  the  system  with  an  effective  energy 
dissipating  mechanism  (radiation  damping)  for  this  particular 
earthquake  accelerogram. 

Similar  conclusions  are  derived  for  the  overturning  moment 


TIME  [sec] 


Fig.  4.  Input  ground  motion  accelerogram:  Taft  earthquake,  N21E 
component  (July  1952),  as  recorded  and  base-corrected  in  the  interval 
f&[6. 014, 7.425]  s.  Peak  ground  acceleration  (PGA)  =  0. 1601  g 
(Input  accelerogram  segment  was  extended  with  zero  values  for  a 
total  duration  of  2.0  s). 


response  at  the  base  of  the  dam,  as  shown  in  Fig.  6,  where 
the  nondimensional  function  H0,  given  by  Eq.  (AAb),  has  been 
plotted. 


Response  Spectra 

The  maximum  absolute  values  of  the  results  provided  by  Eqs. 
(Alb)  and  (44 b)  during  the  history  of  response  for  a  specific  input 
earthquake  may  be  used  to  construct  a  family  of  seismic  response 
spectra  for  base  shear  force  and  overturning  moment,  respec¬ 
tively,  useful  for  dam  analysis  and  design.  The  independent  ab¬ 
scissa  in  the  plot  is  the  reservoir  height,  and  the  parameter 
identifying  each  spectrum  in  the  family  is  the  impedance  ratio  of 
the  sub-bottom  materials  with  respect  to  water.  The  ordinates  rep¬ 
resent  the  extreme  values  of  d0  and  H0  for  the  normalized  base 
shear  force  and  overturning  moment,  respectively.  For  example, 
Figs.  7  and  8  compare  these  spectra  for  a  rigid  dam  when  sub¬ 
jected  to  the  selected  portion  of  the  Taft  earthquake  accelerogram, 
N21E  component  (July  1952),  for  a  very  stiff  reservoir  bottom 
(a=  le+6)  and  for  a  highly  gaseous  reservoir  bottom  (a  =  0.1).  As 
expected,  the  hydrodynamic  forces  are  much  larger  for  the  reser¬ 
voir  foundation  with  a =0.1  than  for  the  rigid  sub-bottom  case. 

Figs.  9  and  10  show  the  comparison  for  the  same  parameters 
of  response,  but  this  time  between  the  spectra  for  the  case  of  a 
rigid  sub-bottom  and  the  case  for  an  “extended  water”  at  the 
reservoir  base.  The  results  suggest  a  substantial  attenuation  of  the 
earthquake-induced  hydrodynamic  forces  (by  a  factor  of  1/4  to 
1/2)  for  reservoir  heights  larger  than  about  45  m  and  amplification 
of  the  base  shear  force  for  reservoir  heights  smaller  than  about  30 
m,  for  the  specific  selected  earthquake,  and  when  the  elasticity  of 
the  solid  sub-bottom  materials  is  included  in  the  model,  as  com¬ 
pared  to  those  forces  for  a  rigid  reservoir  foundation. 

Conclusions 

In  this  investigation,  a  two-dimensional  time-domain  closed-form 
mathematical  model  for  the  hydrodynamic  forces  on  the  upstream 
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Fig.  5.  (a)  Base  shear-force  response  of  a  rigid  dam  facing  a  head¬ 
water  height  of  h  =  65  m,  when  subjected  to  the  selected  portion  of  the 
Taft  earthquake  accelerogram,  N21E  component  (July  1952),  for  sev¬ 
eral  values  of  the  reservoir  bottom  impedance  ratio  a;  (b)  contrasting 
two  of  the  base  shear-force  responses  in  (a):  for  a  very  stiff  reservoir 
bottom  (a=le  +  6)  versus  the  same  response  for  an  almost- water  res¬ 
ervoir  bottom  (a  =1.0) 


vertical  face  of  a  given  rigid  dam  under  a  specified  horizontal 
ground  motion  accelerogram  was  developed.  The  model  includes 
the  absorption  of  energy  at  the  elastic  reservoir  bottom,  which  is 
characterized  by  the  impedance  ratio  of  the  sub-bottom  materials 
with  respect  to  water. 

The  assumption  of  a  rigid,  i.e.,  nonabsorptive,  reservoir 
boundary  may  lead  to  an  overestimation  of  the  earthquake  re¬ 
sponse  of  the  dam-reservoir-foundation  system.  In  particular,  the 
effective  lateral  load  on  the  upstream  face  of  an  existing  concrete 
gravity  dam  may  be  overestimated  and  the  associated  recom¬ 
mended  retrofit  expenditures  unnecessary.  This  was  proved  true  in 
this  investigation  for  reservoir  heights  exceeding  about  45  m 
when  subjected  to  the  selected  input  earthquake  accelerogram.  By 
contrast,  there  was  base  shear  force  amplification  for  reservoir 
heights  smaller  than  about  30  m,  when  the  elastic  sub-bottom  was 


Fig.  6.  Overturning  moment  at  the  base  of  a  rigid  dam  facing  a 
headwater  height  of  h  =  65  m  when  subjected  to  the  selected  portion 
of  the  Taft  earthquake  accelerogram,  N21E  component  (July  1952), 
for  several  values  of  the  reservoir  bottom  impedance  ratio  a 


included  in  the  model,  as  compared  with  the  results  for  rigid 
sub-bottom. 

On  the  other  hand,  a  reservoir  with  sub-bottom  materials  with 
less-than-water  impedance  is  expected  to  exert  earthquake- 
induced  forces  on  a  rigid  barrier  larger  than  those  in  the  case  of  a 
rigid  reservoir  sub-bottom.  This  force  amplification  grows  with 
decreasing  values  of  the  impedance  ratio  of  the  sub-bottom  ma¬ 
terials  with  respect  to  water,  a,  in  the  range  between  0  and  1. 

The  assessment  of  the  impedance  ratio  of  the  sub-bottom  ma¬ 
terials  relative  to  water  (a)  remains  a  difficult  task.  For  one,  this 
parameter  varies  spatially  across  the  bottom  of  the  impounded 
water,  by  virtue  of  the  irregularity  of  the  foundation  layers.  And 


Fig.  7.  Response  spectra  for  the  shear  force  at  the  base  of  a  rigid 
dam  when  subjected  to  the  selected  portion  of  the  Taft  earthquake 
accelerogram,  N21E  component  (July  1952),  for  a  very  stiff  reser¬ 
voir  bottom  (a=le  +  6)  and  for  a  ‘highly  gaseous’  reservoir  bottom 
(a  =  0.1) 
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Fig.  8.  Response  spectra  for  the  overturning  moment  at  the  base  of  a 
rigid  dam  when  subjected  to  the  selected  portion  of  the  Taft  earth¬ 
quake  accelerogram,  N21E  component  (July  1952),  for  a  very  stiff 
reservoir  bottom  (a=le  +  6)  and  for  a  highly  gaseous  reservoir  bot¬ 
tom  (a  =  0.1) 


RESERVOIR  HEIGHT  [m] 

Fig.  10.  Response  spectra  for  the  overturning  moment  at  the  base  of 
a  rigid  dam  when  subjected  to  the  selected  portion  of  the  Taft  earth¬ 
quake  accelerogram,  N21E  component  (July  1952),  for  a  very  stiff 
reservoir  bottom  (a=  le  +  6)  and  for  an  almost  water  reservoir  bottom 
(a  =1.0) 


for  another,  the  probing  mechanism  must  be  able  to  penetrate  the 
bottom  to  a  sufficient  depth  as  to  provide  samples  representative 
of  the  various  substrata. 

During  the  last  two  decades,  several  investigators  have  con¬ 
ducted  projects  in  the  field  attempting  to  quantify  the  parameter  a 
in  real  reservoir-dam  situations.  The  underlying  experiments  were 
based  on  seismic  analyses  of  reflected  and  refracted  pressure- 
wave  signals  generated  with  small  explosives  (Ghanaat  and  Red- 
path  1995)  or  on  the  scattering  of  sonar-acoustics  emissions 
(McGee  et  al.  1995).  The  objective  of  these  waterborne  seismic 
scattering  surveys  was  to  estimate  an  average  value  of  a  for  a 
dam  location  and  the  associated  dispersion,  represented  by  a  co¬ 
efficient  of  variation.  In  a  project  targeting  seven  specific  concrete 
dams  in  the  U.S.  (Ghanaat  and  Redpath  1995),  the  average  value 
of  a  was  inferred  in  the  range  of  0. 3-5.0,  with  a  coefficient  of 


Fig.  9.  Response  spectra  for  the  shear  force  at  the  base  of  a  rigid  dam 
when  subjected  to  the  selected  portion  of  the  Taft  earthquake  accel¬ 
erogram,  N21E  component  (July  1952),  for  a  very  stiff  reservoir  bot¬ 
tom  (a=le  +  6)  and  for  an  ‘almost-water’  reservoir  bottom  (a=1.0) 


variation  of  about  0.25,  with  results  strongly  dependent  on  the 
actual  amount  of  sediments  deposited  in  the  area  close  to  the 
upstream  face  of  the  nonoverflow  monoliths  of  the  dams. 

Notice  that,  generally,  the  effective  impedance  ratio  of  the  sub¬ 
bottom  materials  is  sensitive  to  the  specific  type  of  input  earth¬ 
quake.  Upon  earthquake-induced  oscillations  with  reversal  of 
loading,  the  initial  stiffness  of  the  sub-bottom  materials  may  de¬ 
teriorate,  thereby  modifying  the  effective  elasticity  of  the  under¬ 
lying  strata  and  their  associated  a-value.  Much  more  research  is 
needed  to  evaluate  the  variation  of  response  spectra  for  design 
parameters  of  response  with  the  characteristics  of  the  input 
ground  motion  for  a  given  reservoir-dam-foundation  site. 
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